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The scaling behaviour of randomly branched polymers in a good solvent is studied in two to 
nine dimensions, using as microscopic models lattice animals and lattice trees on simple hypercubic 
lattices. As a stochastic sampling method we use a biased sequential sampling algorithm with re- 
sampling, similar to the pruned-enriched Rosenbluth method (PERM) used extensively for linear 
polymers. Essentially we start simulating percolation clusters (either site or bond), re- weigh them 
according to the animal (tree) ensemble, and prune or branch the further growth according to a 
heuristic fitness function. In contrast to previous applications of PERM, this fitness function is not 
the weight with which the actual configuration would contribute to the partition sum, but is closely 
related to it. We obtain high statistics of animals with up to several thousand sites in all dimension 
2 < d < 9. In addition to the partition sum (number of different animals) we estimate gyration 
radii and numbers of perimeter sites. In all dimensions we verify the Parisi-Sourlas prediction, and 
we verify all exactly known critical exponents in dimensions 2, 3, 4, and > 8. In addition, we 
present the hitherto most precise estimates for growth constants in d > 3. For clusters with one 
site attached to an attractive surface, we verify for d > 3 the superuniversality of the cross-over 
exponent (p at the adsorption transition predicted by Janssen and Lyssy, but not for d = 2. There, 
we find (j) = 0.480(4) instead of the conjectured = 1/2. Finally, we discuss the collapse of animals 
and trees, arguing that our present version of the algorithm is also efficient for some of the models 
studied in this context, but showing that it is not very efficient for the 'classical' model for collapsing 
animals. 
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I. INTRODUCTION 

Lattice animals (or polyominoes, as they are sometimes 
called in mathematics jlj]) are just clusters of connected 
sites on a regular lattice. Such clusters play an impor- 
tant role in many models of statistical physics, as e.g. 
percolation j^], the Ising model (Fortuin-Kastleyn clus- 
ters, Swendsen-Wang algorithm and even lattice 
gauge theories @ . The basic combinatorial problem as- 
sociated to them is to count the number of different 
animals of N sites. Two animals are considered as iden- 
tical, if they differ just by a translation (i.e., we deal with 
fixed animals in the notation of j^), but are considered 
as different, if a rotation or reflection is needed to make 
them coincide. Thus there are e.g. d animals of iV = 2 
sites on a simple hypercubic lattice of dimension d, and 
d{2d — 1) animals with = 3. 

The animal problem can be turned into a statistical 
problem by giving a statistical weight to every cluster. 
In contrast to percolation, where different shapes acquire 
different weights, all clusters with the same number N 
of sites are given the same weight. This is similar to 
self avoiding walks (SAW). A SAW on a lattice is a con- 
nected cluster of N sites with equal weight on all clusters, 
but with a restriction on its shape: each SAW has to be 
topologically linear, i.e. each site is connected by bonds 
to at most 2 neighbours. No such constraint holds for an- 
imals, thus animals are the natural model for randomly 
branched polymers [Tj. 

In addition to animals (or site animals, to be more 
precise) one can also consider bond animals and lattice 
trees. A bond animal is a cluster where bonds can be es- 
tablished between neighbouring sites (just as in SAWs), 
and connectivity is defined via these bonds: if there is no 



path between any two sites consisting entirely of estab- 
lished bonds, these sites are considered as not connected, 
even if they are nearest neighbours. Different configura- 
tions of bonds are considered as different clusters, and 
clusters with the same number of bonds (irrespective of 
their number of sites) have the same weight Q • Weakly 
embeddable trees are bond animals with tree topology, i.e. 
the set of weakly embeddable trees is a subset of bond 
animals, each with the same statistical weight. Strongly 
embeddable trees are, in contrast, the subset of site an- 
imals with tree-like structure. All these definitions are 
illustrated in Fig. ^ 

Like many other statistical models, animals are char- 
acterized by scaling laws in the limit of large N. It is 
believed that all the above statistics (site and bond an- 
imals, weakly and strongly embeddable trees) are in the 
same universality class (same exponents, same scaling 
functions) which is that of randomly branched polymers. 
The number of animals (i.e. the microcanonical partition 
sum) should scale as ^ 



and the gyration radius as 



(1) 



(2) 



Here ^ is the growth constant or inverse critical fugacity, 
and is not universal. In contrast, the Flory exponent v 
and the entropic exponent 9 should be universal. 

In spite of the obvious similarity to the SAW and 
percolation problems, there are a number of features in 
which the animal problem is unusual: 



• The upper critical dimension is d 
1/4 and 9 ^ 5/2 



8. There, v 
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FIG. 1: (a): A site animal with 8 sites; (b) A site tree 
("strongly embeddable tree"); (c) A bond animal which is 
not a tree; (d) A bond tree ("weakly embeddable tree"). 



• The model is not conformally invariant |lfll |. and 
thus the Flory exponent v is not known exactly in 
d = 2. 

• Using supersymmetry, it has been argued by Parisi 
and Sourlas that the animal problem in d di- 
mensions is related to the Yang-Lee problem (Ising 
model in an imaginary external field) in D = d — 2 
dimensions. Based on this relationship (which is 
now proven rigorously ^2|> using a mapping onto 
the hypercubes problem at negative fugacity [T^ ) 
they argued that Q and v should not be indepen- 
dent, but 

e^{d~2)v + \. (3) 

This implies in particular that = 1 in = 2. In 
addition, they showed that = 1/2 in d = 3. 

• Assuming universality so that scaling of the hard 
squares model at negative fugacity can be inferred 
from Baxter's solution of the hard hexagon model, 
and mapping the hard squares model onto lattice 
animals in 4 dimensions, Dhar obtained Q = 
11/6 for d = 4. Thus one knows the exact values of 
V for d = 1, 3, 4, and 8, but not for d = 2, 5, 6, and 
7. 

• In a series of papers, Janssen and Lyssy [TsllT^lr^ 
studied animals attached to an adsorbing plane sur- 
face. For weak adsorption (high temperature) the 
animals have basically the same structure as in the 
bulk, and the partition sum has the same sca ling , 
Z'j^ ~ ju'^iV-® with /x' /X and 6*' = 61 [Hill Ha. 
For strong adsorption (low temperature) there is 



an adsorbed phase. Janssen and Lyssy argued that 
the cross-over exponent between these two phases 
should be super- universal, (p — 1/2 for all dimen- 
sions d> 2. 

In the present paper we address all these points by 
means of a novel Monte Carlo algorithm which follows es- 
sentially the strategy used in the pruned- enriched Rosen- 
hluth method (PERM) [2l|. This is a recursively (depth 
first) implemented sequential sampling method with im- 
portance sampling (bias) and re-sampling ("population 
control"). It seems that PERM in the present imple- 
mentation is much more efficient than previous sam- 
pling methods for animals and trees. Indeed we shall 
present numerous new estimates for critical exponents 
and growth constants which had previously been mea- 
sured only with much larger error bars or not at all. 

All this holds for athermal animals and trees, i.e. when 
there are no attractive forces between monomers. When 
such forces become strong, a number of different collapse 
phase transitions are claimed to occur, depending on the 
detail of the model ^,_2a,„24, 25, 26, 27, 28, 29, 30, 3^ 
HE IMl • The simplest one of these involves site ani- 
mals and a simple contact en ergy for each occupied near- 
est neighbour pair 31, 33. 133 . 13^ and is undisputed. But 
another transition, between two collapsed phases with 
different densities of bonds IHIHIll IH [s^l, is still 
controversial [2^ [2^ |23| . At present all versions of our 
algorithm become inefficient for the first model, when the 
collapse point is approached. This is a bit disappointing 
in view of the fact that PERM for linear polymers is 
dramatically more efficient at the coil-globule transition 
than for athermal SAWs 21] . Obviously this leaves much 
room for further improvements. On the other hand, our 
method should work well for the other transition in large 
parts of the phase diagram. 

Details of the algorithm for site animals will be given 
in Sec. 2. Detailed studies of site animals in the bulk 
and in contact with a wall will be presented in Sees. 3 
and 4. Bond animals and trees will be discussed in Sec. 5. 
Finally, in Sec. 6 we will study animal (and tree) collapse 
due to attractive forces between monomers. The paper 
ends with conclusions in Sec. 7. 



II. NUMERICAL METHODS 
A. Previous Methods 

1. e-Expansions 

Field theoretic e-expansions (where e is the distance 
from the critical dimension) were applied already very 
early to animals and to the Yang-Lee problem [ssf . 
When the relationship between both problems was es- 
tablished, the latter gave the most precise predictions for 
critical animal behaviour in high dimensions {d> 5). 
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2. Exact Enumeration 

Exact enumeration of animals and trees is surprisingly 
non-trivial [s^ [s^, IH, H^- Nevertheless, very exten- 
sive enumerations have recently been performed by I. 
Jensen |EE3j0| for site animals and site trees in d = 2. 
At present they give the best numerical verification of 
the prediction = 1, and they give the most precise es- 
timates for the Flory exponent [y — 0.64115 ± 0.00005) 
and for the growth constants: [i = 4.0625696 ±0.0000005 
for animals and 3.795254 ± 0.000008 for trees. These 
values are more precise than old estimates obtained by 
finite size scaling using strip geometry |42j |. There are 
also enumerations of various animals and trees in higher 
dimensions @, El El H El El 113 E El , but they are 
much less complete and in general they do not at present 
give the best estimates for critical parameters. 

3. Markov Chain Monte Carlo Methods 

The latter are obtained nowadays by Markov chain 
Monte Carlo (MCMC) methods. Such algorithms have 
been used for animals since at least 20 years fs^ , Isol 
m HI]. At present, the most efficient MCMC algo- 
rithm for lattice trees is a version of the pivot algorithm 
HEllHillsll. These simulations showed that = 1/2 
in d = 3, as predicted by Simulations of animals at- 
tached to an attractive wall verified that indeed 0=1/2 
in c? = 2 (as also verified with transfer matrix and 
similar methods [s^ Issjl'). although simulations in d = 3 
gave (j) « 0.714 [5l|, in gross violation of the Janssen- 
Lyssy prediction. 

When applied to SAWs, the pivot algorithm works by 
choosing a pivot point and proposing a rotation of the 
shorter arm around the pivot, and accepting it when this 
leads to no violation of self avoidance ;60] . When adapted 
to trees, one again chooses a random pivot point, but now 
the entire branch hinging on this pivot is cut and glued 
somewhere else. Again this move is accepted only if this 
leads to no violation of self avoidance and if it would not 
lead to wrong cluster topology. 

This method also allows to estimate growth constants, 
if it is used together with the atmosphere method [55l |. 
In the latter, it is counted how many possible ways there 
exist to grow the cluster by one further step, giving in this 
way an estimate o f Zi \[+-\ /Zm- Basically the same method 
had been used in |63| to obtain very precise estimates for 
the critical percolation thresholds in high dimensional 
lattices. 



4- Cluster Growth ('Sequential Sampling', 'Static') Methods 

The first stochastic growth algorithm for trees seems 
to have been devised by Redner [63 ■ Similar methods 
were then used by Meirovitch [6l| and Lam E3| for an- 
imals. But already Leath seems to have realized that 



his well known algorithm for growing percolation clusters 
[6^ could be used also for the study of animals, simply 
by reweighing the clusters. Recently this was taken up 
systematically in |65l |. 

In the following we shall discuss the latter in some de- 
tail, and we shall restrict our discussion to site animals. 
The authors of 65] basically use a standard growth al- 
gorithm for percolation clusters [6^ 66, 67], and then 
reweigh the cluster so that they obtain the correct weight 
for the animal ensemble. In a percolation cluster growth 
algorithm for site percolation, one starts with a single 
seed site and writes it into an otherwise empty list of 
'growth sites'. Then one recursively picks one item in 
the list of growth sites, removes it from the list and adds 
it with probability p to the cluster, and adds all its wet- 
table neighbours to the list. This is repeated until either 
the cluster size exceeds some fixed limit (in which case 
the cluster is discarded), or until the list is empty. A 
cluster with TV sites, b boundary sites, and with fixed 
shape is obviously obtained with probability 

Pm^P^'il-p)', (4) 

i.e. with the correct probability so that any unweighed 
average is just an average over the percolation ensemble. 
Repeating this many times, the animal partition sum is 
then 

^A. = (1/PW=P-^((1-P)-'') . (5) 

The authors of 65] called their method a Rosen- 
bluth method, in view of the obvious analogy with the 
Roscnbluth-Roscnbluth method 68] for SAWs. In the 
latter one also samples from a biased ensemble and then 
reweighs each configuration with the inverse of the sam- 
pling probability to obtain the correct partition sum. 

B. PERM 

Like the original Rosenbluth-Rosenbluth method for 
SAWs, the method of (i^ has the disadvantage that the 
weights have a very wide distribution for large TV. Thus 
even a very large sample will finally, when TV gets too big, 
be dominated by a single configuration, and the method 
becomes inefficient even though it is easy to generate 
huge samples. 

PERM (or any other strategy with resampling) tries to 
avoid this by trimming the width of the distribution of 
weights, by pruning configurations with very low weight 
and making clones of high weight configurations which 
then share the weight among themselves. In many situa- 
tions this has proven to be extremely efficient |6a[70l[7ll ]. 
But we cannot yet apply it to animals, since we have to 
be able to estimate the weight of a cluster while it is still 
growing, and up to now we have only discussed the rela- 
tionship between animals and percolation clusters after 
they had stopped growing. 

In the following we shall again discuss only site ani- 
mals, bond animals and trees being discussed in Sec. 5. 
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FIG. 2: Growing clusters generated breadth first (top) and 
depth first (bottom). In both cases we used p = Pc = 0.5927, 
and in both cases TV = 4000. Occupied sites are depicted by 
smaU points, growth sites by heavy dots. Both figures are 
plotted with the same scale. 



To obtain the relationship between still growing perco- 
lation clusters and animals, let us consider a cluster with 
N sites, g growth sites, and b sites which definitely belong 
to the boundary. At each of the growth sites the cluster 
can grow further, or it can stop growing with probabil- 
ity 1 — p. Thus this cluster will contribute with weight 
(1— p)f to the sample of percolation clusters with N sites 
and b + g boundary sites. Its contribution to the animal 
ensemble is smaller by a factor [p^ {I — p)''^^]^^ , and we 
have thus 

ZA.=p-^((l-p)-^) . (6) 

This is exactly the same formula as above, but now the 
average is taken over all growing clusters, while before we 
had averaged over clusters which had stopped growing. 

Before we can implement these ideas, we have to dis- 
cuss two problems: 

(a) How are the clusters to be grown precisely? 

(b) On what basis should we decide when to prune and 
when to clone a cluster? 

As we shall see, both questions are not trivial. 

(a) Percolation cluster growth algorithms can be depth 



first or breadth first. In the former, growth sites are writ- 
ten into a first-in last-out list (a stack). A typical code 
for this is given in |67| . For a breadth first implementa- 
tion we use, instead, a first-in first-out list ('fifo-list' or 
queue) . Two 2-d clusters growing according to these two 
schemes, withp = Pc = 0.5927, are shown in Fig.|21 Both 
have TV = 4000. But the cluster grown using a stack has 
a completely different shape and has ~ 3 times as many 
growth sites as the one grown with a queue! Most of these 
growth sites are nearly dead (their descendents will die 
after a few generations), but this is not realized because 
they are never tested. Since also the fluctuations in the 
number of growth sites are much bigger in a depth first 
implementation, the weights in Eq.@ will also fluctuate 
much more, and we expect much worse behaviour. This 
is indeed what we found numerically: Results obtained 
when using a stack for the growth sites were dramatically 
much worse than results obtained with a queue. 

Notice that this is independent of the way how pruning 
and cloning is done. Indeed we implemented this "pop- 
ulation control" recursively as a depth first strategy, as 
was done for all previous applications of PERM. 

In addition, there are also some minor ambiguities in 
percolation cluster growth algorithms, such as the order 
in which one searches the neighbours of a growth site 
and writes them into the list. In 2 dimensions one can 
e.g. use the preferences east-south-west-north, or east- 
west-north-south, or a different random sequence at every 
point. We found no big differences in efficiency. 

(b) In most previous applications of PERM, the best 
strategy was to base the decision whether to prune or 
branch directly on the weight with which the config- 
uration contributes to the partition sum |72|. This 
would mean in the present case that we clone, if Wn = 
p~"(l — p)~^ > c+Zn where c+ is a constant of the order 
1 ... 10 and Zn is the current estimate of Z„. Similarly, 
a cluster would be killed (with probability 1/2 |^), if 
Wn < C-Zn with c_ slightly smaller than 1. 

In the present case this would not be optimal, since it 
would mean that mostly clusters with few growth sites 
are preferred (they tend to have larger values of b, for the 
same n), and these clusters would die soon and would 
contribute little to the growth of much larger clusters. 
Thus we defined a fitness function 

/„ = iy„/(l-p)"f =p-"(l-p)-''-"9 (7) 

with a parameter a to be determined empirically, and 
used 

/n>C+(/„), /„<C_(/„) (8) 

as criteria for cloning and pruning. We checked in quite 
extensive simulations that best results were obtained 
with a — 1 (except when TV is small), and in the fol- 
lowing we shall use only this choice. 

Finally we have to discuss the optimal values of p. It 
is clear that we should not use p > pc, where Pc is the 
critical percolation threshold. Since minimal reweighing 
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FIG. 3: Statistical errors of in Z for lattice animals in d = 2 for 
various values of p. To make the different runs comparable, 
errors are multiplied by the square root of the CPU time 
measured in seconds. In 2 dimensions pc = 0.5927. 
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FIG. 4: Same as Fig. 01 but for d = 8 where Pc = 0.0752. The 

straightest curve corresponds to the largest p. 



is needed for small p (subcritical percolation is in the 
animal universality class), one might expect p ^ Pc to be 
optimal. This is indeed true for small values of N (which 
we are not primarily interested in), but not for large N . 
For the latter it is more important that clusters grown 
with p Pc have to be cloned excessively, since they 
otherwise would die rapidly in view of their few growth 
sites. 

To decide this problem empirically, we show in Fig. |3 
the errors of the estimated free energies — — 
for d = 2. More precisely, we show there one standard 
deviation multiplied by the square root of the CPU time 
(measured in seconds), for different values of p. Each 
simulation was done on a Pentium with 3 GHz using the 
gcc compiler under Linux, and each simulation was done 
for iVmax = 4000 (although we plotted some curves only 
up to smaller N, omitting data which might not have 
been converged). We see clearly that small values of p 
are good only for small N. As N increases, the best re- 
sults were obtained for p ^ pc- The same behaviour was 
observed also in all other dimensions, and also for ani- 
mals on the bcc and fee lattices in 3 dimensions (data not 
shown). As an example we show in Fig. ^ the analogous 
results for d = 8. There we used a 64 bit machine (a 600 
MHz Compaq ALPHA), because this simplified hashing 
(for large d we used hashing as described e.g. in 61]). 

Notice that the errors shown in Fig.^jarc much smaller 
than those in Fig.|2| although the machine was slower and 
the animals were twice as large (A^max = 8000). Indeed, 
the errors decreased monotonically with d, being largest 
for d = 2. Using p slightly smaller than pc we can obtain 
easily very high statistics samples of animals with several 
thousand sites for dimensions > 2. A typical 2-d animal 
with 12000 sites is shown in Fig. |31 and a 3-d animal on 
the bcc lattice with 16000 sites is shown in Fig. El 

To check the reliability of our error bars we looked at 
distributions of tour weights as described in ItBI . A tour 
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FIG. 5: A typical lattice animal with 12000 sites on the square 
lattice. 



is the set of all configurations generated by cloning from 
one common start and therefore possibly being strongly 
correlated. If the distribution P of tour weights W is 
very broad, we are back to the problem of the Rosen- 
bluth method that averages might be dominated by a 
single tour. To check for this, we plot P{\nW) against 
In W, and compare its right hand tail to the function 
1/W. If the tail decays much faster, we are presumably 
on the safe side, because then the product WP{ln W) has 
its maximum where the distribution is well sampled. If 
not, then the results can still be correct, but we have no 
guarantee for it. 

In Fig. (71 we show these tour weight distributions for 
two-dimensional animals with 4000 sites, for p = 0.57 
and for p = 0.47 [tH. We see that the simulation with 
p = 0.57 is on the safe side, but not that for p = 0.47. 



FIG. 8: Deviations of the logarithm of the number of 2-d site 
animals from exact enumerations of I. Jensen |4l| . 



FIG. 6: A 3-d lattice animal with 16000 sites on the bcc 
lattice. 
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FIG. 7: Log-log plot of distributions of tour weights of 2d an- 
imals with iV = 4000, for p = 0.57 and for p = 0.47, together 
with a straight line representing the function y = const/W. 
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FIG. 9: Relative deviations of the squared radii of gyration 
from exact values of [4H. 



different ways, plotting effective exponents against dif- 
ferent powers of trying different ansatzes for higher 
order correction to scaling terms, etc. They are not based 
on simply making least square fits over fixed intervals of 
N, as this could lead to very large underestimations of 
corrections to scaling. All quoted numbers are such that 
we believe, to the best of our knowledge, that the true 
value is most likely within one quoted error bar. 

The total CPU time spent on the simulations reported 
in this paper is « 25, 000 hours on fast PCs and work 
stations. 



Similar plots for other simulations described in this paper 
showed that all results reported below are converged and 
reliable. 

Error bars quoted in the following on raw data (par- 
tition sums, gyration radii, and average numbers of 
perimeter sites or bonds) are straightforwardy obtained 
single standard deviations. Their estimate is easy since 
clusters generated in different tours are independent, and 
therefore errors can be obtained from the fluctuations 
of the contributions of entire tours (notice that clusters 
within one tour are not independent, and estimating er- 
rors from their individual values would be wrong). 

On the other hand, errors on critical exponents and on 
growth constants are obtained by extrapolation. This is 
an ill-posed problem, and therefore any error obtained 
this way is to some degree subjective. All such errors 
quoted in the following are based on plotting the data in 



III. SITE ANIMALS IN 2 TO 9 DIMENSIONS 
A. d = 2 

Before we report our final results, we show one more 
test where we compare our raw data for d = 2 with the 
exact enumerations of In Fig. |H1 we show the true 
relative errors of our estimates of the partition sum. Al- 
though there is some systematic trend visible, this is still 
within two standard deviations and thus not significant 
(notice that our values for different N are not indepen- 
dent). Relative errors of the squared gyration radii are 
shown in Fig.|5| These data show that our estimates are 
basically correct, including the error bars. 

Plotting directly our values of Zjv would not be very 
informative, neither would be a plot oflnZN~aN, where 
a = In /i. Both ways of plotting would hide any statistical 
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FIG. 10: Plot of IuZn -aN + einN against In for d = 2, 
with a = 1.401812, 1.4018155, and 1.401819 (top to bottom). 
Error bars are plotted only for the central curve. 
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FIG. 12: Same data as in Fig.Cni but plotted against 1/N° 
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FIG. 11: Plot of R^T^/N"^" against IniV for v = 0.6406, 0.6412, 
and 0.6418. Again, error bars are plotted only for the central 
curve. 



errors. A more meaningful way of plotting our full data 
for Zm is used in Fig. 1101 where we plot \nZM — aN+\aN 
against N for three values of a. Error bars are shown only 
for the central curve, although all three curves have of 
course the same errors. In view of Eq. J^l, and accepting 
the prediction that 9 — 1 , we would expect a curve which 
becomes horizontal for large TV. This is indeed seen for 
the central curve, but the obvious corrections to scaling 
make a precise estimate of /i difficult. The same is true for 
the gyration radii. In Fig.^]we show i?^ /N"^^ against N 
for three candidate values of v. Again strong corrections 
to scaling are seen. 

For these corrections one expects 



and 



~ ii'^N-^l + bzN-^ + ■■■) (9) 



Rl^ N^^il + hRN-^ + ■■■) , (10) 
where A is the correction exponent 9|, hz and are 
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FIG. 13: Same data as in Fig.^l but plotted against l/iV° 



non-universal amplitudes, and the dots stand for higher 
order terms in 1/A^. Notice that A is universal, and is the 
same in both equations. There are several methods dis- 
cussed in the literature for estimating A. We estimated 
it by plotting \iiZm - aN + \yiN and R%/N'^'' against 
X = . Straight lines are expected near x — if and 
only as — A. We could not find a value of S where these 
lines were absolutely straight for all x, but the straight- 
est behaviour near x « was obtained with 6 w 0.9, see 
Figs, inland 1121 

We thus conclude that A — 0.9 ±0.1 which suggests 
that the leading corrections to scaling are analytic (A = 1 
exactly). This is in agreement with the exact enumera- 



tions of |^. l40lf4l| and with the exactly known correction 
exponent for the Lee- Yang problem |3, 123 , but disagrees 
with the Monte Carlo estimate 0.65 ±0.20 of H^. Notice 
that originally ,llj the connection between the Lee- Yang 
and animal problems was established only for the lead- 
ing terms, and therefore the authors of suggested not 
to use the Lee- Yang correction to scaling exponents for 
animals. But the recent proof of |12| | gives an exact map- 
ping between two models in the respective universality 
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TABLE I; Main results for site animals. For convenience we also give in the second column the critical p-values for site 
percolation. 



d 


Pc 


a = In ^ 


e 


V 


e 


A 


Method 


2 


0.5927 


1.4018155(30) 




0.6412(5) 


- 


0.9(1) 


present work 






1.401815696(5) 




0.64115(5) 


- 


1.0 


Series [^,42] 










0.642(10) 


- 


0.65(20) 


MC 


3 


0.3116 


2.1218588(25) 


3/2°).'=) 


1/2").'=' 


]^<2),C) 


0.75(8) 


present work (partially constr.) 






2.1218592(20) 


3/2").'=) 


1/2'"'''=^ 




-j^a),c) 


present work (constrained) 






2.120(2) 










Series [78j 








1.502(3) 








Series [77j 










0.498(10) 




0.54(12) 


MC [56| 


4 


0.1968 


2.587858(6) 


1.835(6) 


0.4163(30) 


1.001(7) 


0.57(8) 


present work (unrestricted) 






2.5878583(40) 


1.833(5) 


0.4181(25) 


0.998(4) 


5/6")'"' 


present work (partially constr.) 






2.5878583(30) 


1.834(4) 


0.417(2) 


■j^a),c) 


5/6")'"' 


present work (partially constr.) 






2.5878483(30) 


11/6").'=) 


5/12").'=) 




5/6").'=) 


present work (constrained) 








1.839(8) 






5/6").'=) 


Series [77] 






2.6012(15) 










MC [42 










0.415(11) 




0.46(11) 


MC [56] 


5 


0.1407 


2.9223194(60) 


2.080(7) 


0.359(4) 


1.001(9) 


0.47(7) 


present work (unrestricted) 






2.9223205(30) 


2.0815(60) 


0.3605(20) 


■j^a),c) 


0.622'')'"' 


present work (constrained) 








2.0877(25) 






0.622(12) 


Series [77] 








2.0807 








e-expansion ^7| 








2.10(3) 


0.367(11) 


■j^a),c) 


0.65(15) 


Series [9] 






2.899(9) 










expansion in l/(2d — 1) [79| 






2.940(15) 










MC [47] 










0.359(11) 




0.40(14) 


MC [56] 


6 


0.1090 


3.1785245(40) 


2.261(12) 


0.315(4) 


1.000(12) 


0.39(6) 


present work (unrestricted) 






3.178521(3) 


2.256(8) 


0.314(2) 




0.412').'=) 


present work (constrained) 








2.2648(15) 






0.412(8) 


Series [77] 








2.2649 








e-expansion 








2.30(4) 


0.325(10) 


■|^a),c) 


0.5(2) 


Series [9] 






3.172(3) 










expansion in l/(2ci — 1) [79| 






3.20(2) 










MC [47] 










0.321(19) 




0.34(13) 


MC [56] 


7 


0.0889 


3.384080(5) 


2.40(2) 


0.282(5) 


0.996(20) 


0.26(6) 


present work (unrestricted) 






3.384079(3) 


2.390(9) 


0.278(2) 




0.205'')'"' 


present work (constrained) 








2.402(5) 






0.205(5) 


Series [^2 








2.4999 








e-expansion 








2.41(3) 


0.282(6) 


■j^a),c) 


0.4(2) 


Series [9] 






3.382(1) 










expansion in l/(2d — 1) [79| 






3.41(1) 










MC [47] 










0.291(11) 




0.35(7) 


MC [56| 


8 


0.0752 


3.554827(4) 


5/2«).<=) 


1/4").'=) 


■j^a),c) 


(+logs) 


present work 






3.5544(7) 










expansion in l/(2d — 1) [7^ 


9 


0.0652 


3.700523(10) 


5/2«).<=) 


1/4°).'=) 


■j^a),c) 


0.25(5) 


present work 



"' Exact value 
*' From III 

"' Used as constraint in the fit 
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1.2 



1.195 



1.19 



1.185 



1.1£ 



0.001 0.002 0.003 0.004 0.005 0.006 

1 /N^ 

FIG. 14: Average number of boundary sites per cluster site for 
2-d animals, plotted against . To reduce finite size effects 
we actually subtracted 4.445 units from (b), before dividing it 
by A'^. Statistical errors are much smaller than symbol sizes. 



classes, and therefore we should use the mapping also for 
the corrections to scaling. 

The critical exponent v and the growth constant can 
be read off Figs. El El and are reported in TablejH The 
latter contains also our main results for all other dimen- 
sions. We see that our estimates for /i and v are still 
much worse than the results obtained by the extremely 
long exact enumerations of Jensen, but they are more 
precise than all other previous estimates. 

We have also made "unbiased" fits where we did not 
assume the theoretical values 9 = 1. We do not show 
details, we just mention that our data would seem to 
exclude 16*- 1| > 0.002. 

Finally, we show in Fig. 1141 the average numbers of 
boundary sites. More precisely, with (6) being this aver- 
age, we plot ((6) - 4.445)/A^ against Subtracting 
4.445 units was done in order to reduce finite size effects. 
Without the very large correction 4.445/iV, this term 
would dominate any other correction term, and would 
mask in particular any possible non-analytic correction. 
The linear shape of the curve suggests that there are no 
non-analytic corrections, and that the next to leading 
term is ~ but the data are too poor to allow a 

firm conclusion. 



B. d > 2 

For d = 3 we show the data for Zjv and in Figs. 1151 
and El plotted in the same way as in Figs. El and El 
Now the straightest line is clearly obtained for J < 1, i.e. 
there definitely seem non-analytic corrections to scaling. 
The best fit was obtained with A = 0.75 ± 0.08 (upper 
panels in Figs.|15|and ll6|l . But from the Lee- Yang prob- 
lem we know 77] that we have A = 1 also in d = 3. 
As seen from the lower panels in Figs. 1151 and 1161 where 
these data are plotted against 1/A^, this is clearly not sup- 



3 = 2.1218563 
a = 2.1218588 
3 = 2.1218613 




0.07 



0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 
1 /N 



FIG. 15: Similar to Fig. [El but for d = 3. The straightest 
curve was now obtained by plotting the data against l/N'^''^^ 
(upper panel). In the lower panel, the same data are plotted 
against 1/N. 



ported by our data. But we cannot, of course, exclude the 
possibility that this is due to very large higher order cor- 
rections to scaling. In view of this we show in TableQltwo 
fits, one unrestricted where A is fitted from the present 
data and one constrained fit where A = 1 is imposed. In 
both fits the Parisi-Sourlas condition 6=1 + {d— 2)^ and 
the exact values i/ = 1/2 and 9 = 3/2 were also used as 
constraints. Fits without imposing v = 1/2 and 9 = 3/2, 
and without assuming the Parisi-Sourlas relation, gave 
bigger errors for the growth constant, but gave exponents 
in full agreement with the predictions: v = 0.500 ±0.002 
and 9 = 1.500 ± 0.001. 

The problem with the correction to scaling exponent 
is obviously due to large sub-leading corrections. It per- 
sists also in higher dimensions. For d = 4, e.g., we 
estimated A = 0.57 ± 0.08, while the exact value ob- 
tained from the Lee- Yang problem is A = 5/6. We 
present therefore in Table ^ four different fits with var- 
ious constraints: One completely unrestricted, another 
with A = 5/6 imposed, a third with the Parisi-Sourlas 
relation imposed in addition, and a final one with even 
the values 9 = 11/6 and — 5/12 fixed. Notice that 
the growth constant can be obtained without knowing 
9, if the Parisi-Sourlas relation is assumed: In this case 
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0.29 



0.28 



0.27 



0.26 



0.25 





V = 0.4985 




v = 0.5 




v = 0.5015 











0.02 



0.04 0.06 

1 / n''-^^ 



0.08 



0.1 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 
1 /N 



FIG. 16: Similar to Fig. 1131 but for d = 3. In the upper panel, 
B?j^/N^" is plotted against l/N^'"^"', in the lower panel it is 
plotted against 1/A'^. 



From the values 



NZkR"!^^ fi^il + const /N"^ + . 
listed in Table HI we see that all four fits are mutually 
consistent. 

Similar fits were also made for d = 5,6 and 7, but we 
list in Tabled only the results of the unrestricted and of 
the completely restricted fits. In all cases the agreement 
between the fits is very good, showing the consistency of 
the data. 

For d = 8 we show in Figs. El and E] the data for 
free energies and for gyration radii, plotted against In N. 
From these plots one sees clearly that the data agree with 
the predicted exponents 6 = 5/2 and v = 1/4. But there 
are very large (presumably logarithmic) corrections, com- 
patible with the fact that d = 8 is the upper critical di- 
mension. We have not tried to make a detailed fit to 
these corrections, since we are not aware of any theoret- 
ical prediction beyond the leading order , and because 
verifying logarithmic corrections is notoriously difficult. 

In Fig. El the average number of boundary sites are 
plotted in a way similar to Fig. 1141 This time an even 
bigger term 6.5/iV had to be subtracted, in order to see 
any possible non-analytic term. The fact that the curve 
is reasonably straight when plotted against l/N"^ sug- 
gests again (as for d = 2) that there is no non-analytic 




10000 



FIG. 17: In Zn - aN + 2.51niV for d = 8, plotted against 
IniV, with three different values of a: 3.554822, 3.5548327, 
and 3.554832 (top to bottom). Error bars are again plotted 
only for the central curve. 




10000 



FIG. 18: Plot of Rf^/VN against InN for animals in d = 8. 
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1 /N^ 



FIG. 19: ({&) — 6.5)/A'^, where (b) is the average number of 
boundary sites per cluster site for 3-d animals, plotted against 
I/N"^ . Statistical errors are much smaller than symbol sizes. 
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(0-1 )/(d-2) 



0.65 
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0.35 



0.25 

23456789 
dimension d 

FIG. 20: The critical exponents v and (6' — — 2) against 
d. 



0.02 
0.015 

0.01 
0.005 


-0.005 
-0.01 



\ \o(1/a^) 






' 0(lV)\ 




- 









0{1/a^) 


0(1/a^) 


... ■■Q 


5 6 


7 8 


9 



FIG. 21: Difference between the measured values of In /.i(ci) 
and successive truncations of the expansion Eq. H12^ , plotted 
against d. Error bars are smaller than the sizes of the symbols. 



TABLE II: Asymptotic ratios between boundary and cluster 
sites, limjv^oo {&)/A'' (column #2); ratios between partition 
sums of bulk and wall grafted animals (column #3); and av- 
erage number of contact of wall-grafted animals with that wall 
(column #4). 



dimension 


limjv-.oo{&)/iV 


limjv-,oo Z'j^/Zn 


(m) 


2 


1.1951(1) 


1.987(8) 


2.892(2) 


3 


2.7877(1) 


2.97(3) 


5.07(1) 


4 


4.5859(2) 


3.98(5) 


7.50(5) 


5 


6.4909(2) 


4.91(4) 


10.12(6) 


6 


8.4503(1) 






7 


10.4363(2) 






8 


12.4346(1) 






9 


14.4378(2) 







correction term. Our estimates for the critical exponents 
and for the growth constant are given in Table |l| 

The estimates for v and for {9 — \)/{d — 2) obtained 
by the unrestricted fits are also shown in Fig. 1201 Ac- 
cording to Parisi and Sourlas, they should coincide. The 
agreement is practically perfect. Moreover, both esti- 
mates have roughly the same errors, and estimating v 
indirectly, using Eq. J^Jl, seems to give slightly smaller 
errors for d> b than the direct measurement. For d — Q 
and d — 7 our results are in very good agreement with 
the e(= 8 — d)-expansion results of • Our results are 
also in very good agreement (< la) with the series ex- 
pansion of the binary Gaussian molecule mixture of |77j | , 
which gives for high dimensions the most precise previous 
exponents for the Lee- Yang problem. 

The growth constants seem to grow linearly with di- 
mension, 



lj{d) w 5.49d- 8.94 



(11) 



for large d, although there are small but statistically 
significant deviations. More precise expressions for the 



large-d behaviour are obtained from an expansion of In fi 
in powers of l/cr with a = 2d — 1 : 



In /i(c?) 



In (7 



2 

1 - - 
a 

79 317 



18321 123307 



24cr2 24cr3 320a^ 



2400-5 



(12) 
0(a-6). 



A comparison of our data with different truncations of 
this expansion is shown in Fig. 1211 This comparison sug- 
gests strongly that the expansion is only asymptotic: For 
any fixed d, there is an order beyond which the expansion 
gives values smaller than the true value, and continues 
to decrease with increasing order. 

Our estimates for the asymptotic number of boundary 
sites per cluster site are given in Table Hll The latter were 
all obtained by assuming no non-analytic corrections to 
scaling, since our data can be fitted for all d to {b)/N = 
(3 + (ii/N + P2/N^' with A' w 2. For large d, our data 
seem to scale as 



(3{d) = lim {b)/N' 

N^OO 



2d — const. 



(13) 



IV. ANIMALS ATTACHED TO A WALL 
A. Athermal Walls 

In this section we will consider rf-dimensional animals 
grafted with one monomer to an impenetrable planar wall 
modeled by a hypcrplane Xd = 0. For this case, it was 
shown in [iq that the partition sum, written now Z'j^ 
instead oi Zn, scales as 



with the same as in the bulk, and with [2C 



(14) 



(15) 
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The last equation has a very simple heuristic explana- 
tion. Let us first map lattice points x = (a;i . . . Xd) on a 
lattice of size L'^ onto integers 



= xi + X2L . . . XdL 



d-l 



(16) 



(we actually used this in our codes to index points 
by a single integer, as this simplifies programming and 
makes memory access faster). Consider now the prob- 
lem of counting all animals restricted to the half space 
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FIG. 22; Ratio Z'j^/Zn for 2-d site animals, plotted against 
(lower curve) and average number of contacts of 
grafted 2-d animals (upper curve). Statistical errors are 
smaller than the sizes of the points. 



Animals attached to an attractive surface 



xq belongs to the animal. On the one hand, this means 
just that we consider animals with fixed positions: for 
each shape we consider only that animal whose "small- 
est" point is Xq. Since we had counted only once all 
cluster shapes related by translations, this means that 
the partition sum obtained now is exactly equal to Z^. 
On the other hand this model is equivalent to the an- 
imal being grafted to an impenetrable wall located at 
Xd ~ Xd,o which is however not quite flat: Xd jumps 
from Xdfl + 1 to Xdfl when any one of the other coordi- 
nates Xj {j = 1, . . .d - 1) goes through Xj^Q. The proof The partition sum now is written as 
of ^3 then just shows that the scaling behaviour is in- 
dependent of these steps, and is the same as for a flat 
surface. In addition, this argument shows that > Zn 
for all N 18]. Indeed, one easily sees that the ratio 
{niN)o = Z'j^/Zn is just the average number of con- 
tacts a free animal in the bulk would have with a flat 
imaginary wall placed just below it. This is not equal to 
the average number of contacts (mAc) of a grafted animal 
with its wall, because the latter animals are counted m 
times if they have m contacts: If we denote by ZM,m the 
number of configurations with m sites in the bottom hy- 
perplane and N sites total (so that Z^ = ^M,m), 
and by ^ the analogous quantity for grafted ani- 



N 



z(;'(g)=^AA.(m)g" 



(18) 



mals, then Z^,/,„ — mZM,m- Therefore, the average 
number of contacts of a grafted animal is given by the 
second moment of m in the bulk ensemble divided by 
the first moment, (mAr) = Y.m^Z'j^i,ml Y.mZ'M,m = 

In Fig. 113 we show the ratio Z'j^/Zn against l/iV°-^, 
together with the average number of contacts (m^r) of 
grafted animals, for d = 2. We see again straight lines, 
showing that 



where A]\[{m) is the number of configurations of lattice 
animals with N site having m sites on the walls, and 
q — e'^^^^ is the Boltzmann factor, e > is the attractive 
energy between the monomer and the wall. 

As q ^ 1, there is no attraction between the monomer 
and the wall, i.e. Z^\l) = Z'^. On the other hand 
it becomes clear that any cluster will collapse onto the 
wall, if q becomes sufficiently large. Therefore we expect 
a phase transition from a grafted but otherwise detached 
to an adsorbed phase, similar to the transition observed 
also for linear polymers. 

Exactly at the transition point q = qc we expect the 
usual scaling laws 



and 



R 



N 



(19) 



(20) 



Z'n/Z. 



N 



{rriN) 



const — 



const 



(17) 



with A' = 0.9 ± 0.05. Similar results were found for 
larger dimensions. The values of A' are close to those of 
A but somewhat larger, and we see no theoretical reason 
why they should be the same. We do not quote numbers 
since they are rather poorly determined. The asymptotic 
values of Z'-^jZj^ and (m^r) are given in Table UTI 



In analogy to critical surface phenomena where this tran- 
sition would correspond to the "special" point [s^l, we 
expect V to be the same Flory exponent as in the bulk, 
while ds should be a new and independent exponent. The 
growth constant /i, although being not universal, should 
be the same as in the bulk. 

Away from the critical point we expect a scaling ansatz 



^^{{q-q,)N% 



(21) 



with the crossover exponent being a second new expo- 
nent. Taking the derivative of XvlZ^j!^ {jq) with respect to 
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q and setting q = qc thereafter, we obtain for the average 
energy 



E^iq,) = (em) ~ TV*. 



(22) 



Taking two derivatives we obtain for the specific heat per 
monomer near (but not exactly at) the critical point 



CNiq) 



1 



NKT^ 



{{(emf)-{emf)^{q-q,)-" (23) 



with 



while 



a = 2 - 1/0, 



CN{qc)-N^*~ 



(24) 



(25) 



In principle, all four scaling laws can be used to locate 
the critical value qc- With conventional (Metropolis type) 
Monte Carlo simulations one cannot use easily Eq. (|19|l . 
since precise estimates of the partition sum are difficult 
to obtain. In this case it is usually Eq. (|23l) which is 
used. With PERM we do have very precise estimates of 
Z^^\q), and therefore we can use Eq. ((T^ . but we shall 
see that it is indeed Eq. which gives - together with 
the two others - the most precise estimate. This is very 
similar as for adsorption of linear polymers [U . 

In the following we shall assume e — 1 without loss of 
generality. In order to compare with previous analyses 
we want to have specific heats for discrete values of N , 
but for a continuous range of q. They are most easily 
obtained from histograms 



'm,m' 



(26) 



Confie 



which are normalized such that z''^\q) = "^^Pim^q). 
Notice that we obtain from the simulations not only the 
shape of the histogram, but also its absolute normaliza- 
tion, which makes it easy to combine two histograms ob- 
tained in runs with different nominal values of q. All 
we have to know are rough values of their relative sta- 
tistical errors. These we can estimate from the num- 
ber of tours which contribute to a particular value of m, 
AP{m;q) cx Although this estimate is not 

very precise, it is fully sufficient to obtain smooth global 
histograms by joining histograms which cover narrow re- 
gions in m. 

Specific heats for 2-d animals with lengths up to = 
1200 are shown in Fig. These data are very similar 
to the results of [ESl, although the latter are for trees. 
According to Janssen and Lyssy jl^ 0| we expect 
a — 0, i.e. the specific heat curves for different N should 
intersect exactly at the critical point. This gives roughly 
qc = 2.27. But a close look at the insert in Fig. 1231 re- 
veals that these intersections slightly shift to larger q as 
N increases. Thus we have considerable corrections to 




FIG. 23: Specific heats CAr(q) per monomer for 2-d site ani- 
mals, plotted against q for various values of A'^. Smooth curves 
show results obtained by histogram reweighing, points indi- 
cate results of single runs. The insert shows the region around 
q ~ 2.27 where the curves intersect. 
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FIG. 24: Z^^\q)/tJ," against iV for various values of q, for 
site animals on the square lattice grafted to a wall. Error 
bars here and in the following figures are comparable to the 
thickness of the lines. 



scaling, preventing us from attributing error bars to this 
estimate. 

Alternatively, we turn towards the partition sum itself. 
In Fig. 1^ we show log-log plots of Z^^(g)//i^, for var- 
ious values of q close to qc- We see the expected power 
law, but determining the critical point from this figure is 
difficult because of the substantial corrections to scaling. 
We thus multiply with an estimated power N^" and plot 
the data against 1/A^°-^. The result is shown in Fig. 1251 
(where we actually plot the logarithm on the j/-axis). No- 
tice that we did not have to make a new estimate of /i; 
rather, we could take the old and very precise estimate. 
The value 9s = 0.87 was chosen so as to give the best 
straight line for a suitably chosen qc. Indeed, from this 
plot we would conclude that qc « 2.278. 

We now turn towards Eq. H22|l . In Fig.|2n|we plotted 
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FIG. 25: Part of the data shown in Fig. 1241 after interpolating 
them to a finer set of 5- values, multiplying them with N^'^"^ , 
and taking their logarithm, plotted against 
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FIG. 27: Part of the data shown in Fig.UHl but multiplied by 
iVi-^-* and plotted as curves with fixed q against 1/N ' . 
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FIG. 26: E/N"-^^ against 1/N°-^ for various values of q, for 
site animals on the square lattice grafted to a wall. 



E/N'^-^^ instead of E/\/N. This was chosen because it 
suggests that qc « 2.278, in agreement with the value ob- 
tained above from Zj^. Assuming <j) = 1/2, in contrast, 
would have given qc ~ 2.283 which would be incompati- 
ble with the data for Z^. 

Finally, we plot N'^~'^'f'CN{q) against l/iV°-^, in order 
to compare with Eq. (|25(l . As seen from Fig. [23 this is 
fully consistent with qc ~ 2.278. If we had taken = 1/2, 
we would again get a too large estimate qc ~ 2.285. 

Summarizing, we obtain as our best estimates: 

= 2.2778 ±0.0008 , = 0.480 ± 0.004 (27) 

together with 

61^ = 0.870 ± 0.009 , A, =0.8 ±0.2. (28) 

The large error of reflects the fact that the best esti- 
mates obtained from the different observables would be 
quite different, suggesting again large non- leading cor- 
rections. But this seems to have little effect on the esti- 
mates of the other quantities. Since we believe that we 



have taken into account all systematic errors, we claim 
that the Janssen-Lyssy conjecture (j) = 1/2 is slightly but 
significantly violated in d = 2. The previous estimates 
(j) = 0.505(15) [13 and = 0.503(3) ^ most hkely suf- 
fer from such systematic errors. On the other hand, our 
result is in agreement with the MC estimate </> = 0.50(3) 
of [5^ . Surprisingly, our estimate for (f) agrees within 
the error bars with the most recent estimate of the cross- 
over exponent for unbranched polymers attached to an 
attractive wall in 3 dimensions, (psAW = 0.484(2) [s^ . 
while 0SAW = 0.5 for unbranched 2-d SAWs 0|. 

Before leaving this problem, we should discuss the gy- 
ration radii. Their behaviour near the critical adsorption 
point is somewhat more complicated. For q > qc the 
gyration radius scales as Rn ^ N. At q < qc we ex- 
pect it to scale as Rn ~ N'^ with the same value of ly as 
in the bulk, as in other surface critical phenomena |80j |. 
The effect of the wall is only seen then in the amplitude 
A = \im]\j^ao{Rff) /N'^" . For q < qc it should be larger 
than for free animals, and for q — qc we should expect 
it to be even larger, Abuik < ^waii < ^c- The reason is 
that the main effect of the wall is to squeeze the animal 
in the direction perpendicular to the wall, which by the 
excluded volume effect makes it more extended in the 
direction parallel to the wall. 

This is a bit analogous to the case of an unbranched 
polymer between two athermal walls 's^'l : If the distance 
D between the walls is decreased, at first the shrinking of 
the perpendicular extension dominates any increase par- 
allel to the plates. However, if D is much smaller than 
the Flory radius, the stretching parallel to the walls dom- 
inates, and Rn increases in comparison to a free polymer 




Our data fFigl28|l indicate that Rn/N'^ is larger than 
for animals in the bulk (where Abuik ~ 0.189), and that 
it increases with q. But at q ~ qc it is not monotonic 
in N: It increases with A'' until A'' « 300, and then de- 
creases sharply. This strange behaviour might have been 
expected from the analogy with unbranched self avoid- 
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TABLE III: Critical Boltzmann factors, crossover exponents, 
and critical exponents at the adsorption transition for site 
animals on simple (liyper-)cubic lattices grafted to a flat at- 
tractive wall. 



dimension 









2 


2.2778(8) 


0.480(4) 


0.870(9) 


3 


1.4747(6) 


0.50(1) 


1.476(7) 


4 


1.2674(6) 


0.50(2) 


1.91(1) 


5 


1.1786(5) 


0.51(3) 


2.18(4) 
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FIG. 28: Ratios {R%)/N'^'^, where Rn is the gyration radius 
of grafted 2-d animals, plotted against l/A''"'*. 



ing walks between two athermal walls. For small N and 
q = Qc the effect of the wall is strong, and the increase of 
the size parallel to the wall dominates. But for N oo 
the effect of the wall becomes increasingly weaker, and 
the stretching along the wall becomes less important. We 
verified that it is indeed the slower increase of the parallel 
component which lets Rn/N" decrease for large N, but 
we found no similar effect in simulations (unpublished) 
of unbranched polymers at the critical adsorption point. 
Thus we have at the moment no good explanation for 
this effect. 

For higher dimensions, the same kind of analysis as in 
d — 2 gave the estimates given in Table IIIII The main 
problem in these analyses is again that the best estimate 
for the leading correction to scaling exponents for the 
different observables did not quite agree with each other. 
This hints at the presence of more than one important 
term in the scaling corrections, and it dominates the error 
estimates. The most remarkable result seen in Table Hill 
is the perfect agreement with the Janssen-Lyssy predic- 
tion (/) = 1/2 in all dimensions > 3. In particular, it 
seems that the strong violation seen in ^9] for = 3 was 
due to an underestimation of finite size effects. Actually, 
Janssen and Lyssy had derived (j) — 1/2 only for d = 3, 4, 
and > 8, but not for d = 2. 

Finally, we show in Fig. l29l the monomer density profile 
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FIG. 29: Monomer density profiles p{z) (with z being the dis- 
tance from the wall) for a 3-d animal grafted to an attracting 
wall. 



p{z) for 3-d animals, where z is the distance from the wall. 
Animal sizes for this figure were N = 4000. All curves 
must of course go to zero for z oo. At temperatures far 
above Tc, i.e. at g <C (?c, the entropic repulsion from the 
wall dominates and p{z) has a maximum at a finite value 
of z. In the adsorbed phase p{z)has its maximum at z — 
and decreases monotonically with z. Notice that the 
transition from non-monotonic to monotonic behaviour 
does not happen exactly at qc, but for q slightly smaller 
than Qc- Presumably this is a finite size effect, and the 
transition would happen at qc for much larger animals. 



TREES AND BOND ANIMALS 



A. Site Trees 

The simplest modification of the codes presented so 
far is needed for simulating site trees. As we pointed out 
in the introduction, site trees are site animals without 
loops. Thus the number of nearest neighbour pairs is 
just A^ — 1 for a tree of N sites. It is easy to count 
the number of occupied nearest neighbour pairs as the 
cluster grows. We have just to prune the growth as soon 
as this number is equal or larger than N. Apart from 
that, pruning and branching is done exactly as before, 
and all comments made in Sec. 2 about the efficiency of 
the algorithm apply also to site trees. 

We made simulations only for 2-d site trees, and only 
with rather modest statistics. Our results were fully in 
agreement with those of Jensen |4l| . In particular we 
obtained fi — 3.79527(4) after three days of CPU time 
on a 3GHz Pentium, to be compared to the estimate 
3.795254(8) obtained in 
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B. Bond Animals and Bond Trees 

In order to simulate bond animals and bond trees, 
one has to grow bond instead of site percolation clus- 
ters. Cluster growth algorithms for bond percolation are 
very similar to those for site percolation and about as 
easy. One just has to remember that in bond percolation 
one often does not distinguish between clusters with the 
same configurations of sites, but with different bond con- 
figurations. However, for animals it is essential to make 
this distinction. 

Let us denote by k the number of non-bonded nearest 
neighbour pairs (often called 'contacts' in this context), 
by b the number of surface bonds, and by m the number 
of established bonds between nearest neighbours. In the 
bond percolation ensemble, a cluster with these 'quantum 
numbers' has a weight (cf. Eq.(21)) 

PNkbrn ^P"'il ■ (29) 

This is slightly more complicated than in the site perco- 
lation case, but one can follow the same strategy when 
using this to simulate (bond) animals. We just have to 
replace the number of perimeter sites in the weight fac- 
tor by b + k, and if we want to simulate trees, we have 
of course to prune all configurations which are not tree- 
like. Growth sites have to be replaced by growth bonds. 
Moreover, it is a bit more natural to consider ensembles 
with fixed to, i.e. with fixed number of established bonds, 
rather than sites 0- 

The heuristics worked out in Sec. 2 remain valid: one 
obtains much better results when the trees are grown 
breadth first instead of depth first; one should use a fit- 
ness function /(C) = (1 — p)~^W{C), where g is now the 
number of growth bonds; and one should simulate at a 
slightly subcritical value of p which approaches Pc as the 
trees to be simulated become larger and larger. For the 
same values of N, the optimal values of p/pc were how- 
ever somewhat smaller. For d — 2 and N — 1000, e.g., 
best results were obtained with p « 0.45 (with pc = 1/2). 

Let us first discuss trees. In two dimensions, we ob- 
tained fi = 5.14276±0.00002. This is compatible with the 
best previous estimate, p. = 5.14339 ± 0.00072 but 
more than an order of magnitude more precise. Our esti- 
mate is based on rather small clusters (iVmax = 500), but 
very high statistics: The error of In Zjv is A In Zn^^qo = 
0.000451. These simulations confirmed also that the lead- 
ing correction to scaling exponent is between 0.9 and 1.0. 
In three dimensions we obtained fi — 10.54646 ±0.00010. 
This case is a bit special, since it has strange looking 
corrections to scaling, as seen from Fig. 1301 Apart from 
Fig. I2H1 this is our clearest and most striking example 
showing that the corrections cannot be described by a 
single power, and illustrates the pitfalls in estimating 
correction to scaling exponents from poor data. In other 
cases, the existence of more than one power in correc- 
tions to scaling have often to be inferred less directly, 
e.g. by comparing different observables (see Sec. 4B) or 
by invoking universality |85| . 
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FIG. 30: Plot of In^jv - aN + 3/2 In against l/iV'^ '^^ for 
bond trees on the simple cubic lattice. 



General animals were simulated only with somewhat 
lower statistics, since we are not aware of any other re- 
cent high statistics simulations. On the square lattice we 
obtained n = 5.20789 ± 0.00004, to be compared with 
the previous series expansions estimate 5.208 ± 0.004 of 
[78l|. On the simple cubic lattice we got /i — 10.61539 ± 
0.00006, to be compared with 10.63 ± 0.05 and 
10.62 ± 0.08 [86|. Mean-square gyration radii were e.g. 
R%^^gg = 347.974 ± 0.036 in d = 2 and iJ^^ggg = 
165.669 ± 0.022 in d = 3, to be compared to the best 
previous MC estimates, 348.32 ± 0.88 and 166.03 ± 0.34 
[3 (the quantity n displayed in the first lines of Tables 
Al to A7 of 8] is not the number of bonds, as stated 
there, but one plus this number). 

Finally we point out that we can also use bond percola- 
tion as a starting point for the simulation of site animals. 
We just have to use the fact that site animals are isomor- 
phic to the subset of bond animals with maximal number 
of bonds for a given configuration of occupied sites. Us- 
ing this we obtained for d = 2 results in agreement with 
those of Sec. 3, but the algorithm was somewhat less ef- 
ficient than that based on site percolation. 



VI. ANIMAL COLLAPSE 

A. Collapse of Site Animals 

In order to describe collapsing animals and/or trees, 
one has to introduce attractive monomer-monomer inter- 
actions. Historically the first model of this type |^ [s^ 
starts from site animals and introduces a contact energy 
for each 'contact', where a contact is a pair of occupied 
nearest neighbour sites (notice that the definition of con- 
tacts used here differs slightly from that used for bond 
animals). Let us denote this energy as — e, and the cor- 
responding Boltzmann factor as g = exp(/3e) > 1. Let us 
furthermore denote by to the number of contacts. The 
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partition sum is then written as 

ZN{q) =J2^N,^nq"' (30) 

rn 

where Cjsf^rn is the number of different clusters with 
N sites and m contacts. If the clusters are embedded 
in some solvent, repulsive monomer-solvent interactions 
need not be included explicitly, since the number s of 
solvent contacts satisfies AfN = s + 2m, where Af is the 
coordination number of the lattice {J\f = 2d on a simple 
hypercubic lattice). 

Simulations with PERM are straight forward. We 
just have to modify the weight factors Wn to W]y_,n = 
Wnc^^™. We found that again, as for the previous case 
q = 1, it was better to simulate clusters breadth first 
than depth first. We also verified that it was advan- 
tageous to include a factor « (1 — p)~^ in the fitness 
functions, just as for athermal clusters. But the results 
were rather disappointing, at least for low dimensions. 
This might seem at first surprising, given the fact that 
PERM works extremely well at the collapse transition of 
linear polymers in 3 dimensions (2l| . But it is easy to 
see the reason for this difference. For linear polymers, 
d = 3 is the upper critical dimension, and 9 polymers 
form essentially random walks with very small logarith- 
mic corrections. Thus starting off with random walks, 
PERM can do with very few resampling steps. There is 
only one pruning or cloning needed for every 2000 sim- 
ple forward steps Collapsing site trees, at least in 
low dimensions, are however very different from site per- 
colation clusters. A convenient observable to see this 
difference is the average number (to) of contacts. For 
both models, (to) is roughly proportional to N. But 
for site percolation on the simple cubic lattice one finds 
(to) /N « 0.15, while the same number for collapsing an- 
imals (at q w 3.22 ^) is « 0.40. Thus there is still a 
huge amount of re-sampling needed, even more than for 
athermal animals where {m)/N « 0.065. 

For this reason we do not even show any data for the 
collapse in d = 2 and d ~ 3. The situation improves 
when d is increased. Results for d = 4 are shown in 
Fig. 1^ We see a rather sharp collapse transition at 
q — Qc — 2.98 ± 0.02. It is hard to give precise estimates 
of any critical exponents from these data. But the Flory 
exponent v seems to be the same as for athermal animals, 
within rather small errors. This would not be easy to 
understand, if it were exactly true. Anyhow, we are not 
aware of any theoretical prediction to compare this with. 
We are not even aware of any prediction of the upper 
critical dimension for this collapse transition. 

B. Collapse of Bond Animals and Bond Trees 

Let us now switch attention to collapse models based 
on bond animals. We have now two different possible 
interactions. Instead of the single parameter e for the in- 
teraction strength in site animals, we can now introduce 
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FIG. 31: Squared gyration radii of collapsing site animals in 
four dimensions divided by A^^", where = 5/12 is the exact 
Flory exponent for non-collapsing animals. 



different parameters ei for bonded and £2 for non-bonded 
neighbour contacts, and define y = e*^^ ,t = e*^^. Notice 
that as before we do not have to introduce also an addi- 
tional interaction with the solvent, because the number 
b of surface bonds is not independent of to, k, and N: 

2m + 2k + b = AfN. (31) 

We then define 

ZN{y,T)^j2^^-"^''^y"'^' ■ (32) 

m.k 

The model discussed in the last subsection is obtained 
by taking the limit t — |2^. In this limit only bond 
configurations with the maximal number of bonds (for 
a given configuration of sites) contribute to the parti- 
tion sum. On the other hand, bond percolation without 
reweighing corresponds, due to Eqs. H29|l and H31|l . to the 
curve 

y=p/{l^p)\ r=l/(l-p), 0<p<l. (33) 
or, explicitly, 

y = T(r - 1) , r > 1 . (34) 

Critical bond percolation {pc — 1/2) corresponds to 
y = T = 2. Simulations should be very easy in the neigh- 
bourhood of this point, but they also should be not too 
difficult in the neighbourhood of the entire bond percola- 
tion line. The reason is simply that along the entire line 
one should not need much resampling. This should be 
enough to obtain precise estimates for large parts of the 
phase diagram, and in particular to clarify the existence 
of two different collapsed phases. Results will be given 
in a forthcoming paper. 
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VII. CONCLUSION 

We have shown that the basic idea of PERM, namely 
the recursive implementation of biased sequential sam- 
pling with re-sampling, can be applied also to lattice an- 
imals and lattice trees. These are two models for ran- 
domly branched polymers. Our algorithm is extremely 
efficient (except for the collapse of site animals), obvi- 
ously much more efficient than previous Monte Carlo al- 
gorithms. We applied it to simple (hyper-)cubic lattices 
in up to nine dimensions, but made also less complete 
simulations on bcc and fee lattices. Our algorithm works 
indeed better in higher dimensions, nevertheless we ob- 
tained high statistics results also for large animals in two 
dimensions. 

We verified a number of theoretical predictions. In par- 
ticular, we verified the Parisi-Sourlas connection between 
entropic and Flory exponents, and we verified the values 
of these exponents whenever they are exactly known. We 
also verified that the cross-over exponents for branched 
polymer adsorption on plane walls is super-universal for 
d > 2, as predicted some time ago (but not for d — 2!), 
and we gave precise estimates of the other critical expo- 
nents at this adsorption transition. 



There are a number of problems we did not yet study, 
although our algorithm seems ideally suited for them, 
and which we plan to address in forthcoming papers. One 
is that of collapsing animals where we hope to be able to 
verify or disprove the existence of two different collapsed 
phases. Another is the dependence on the wedge angle, 
of the entropic exponent of a 2-d animal grafted at the tip 
of this wedge. In conformally invariant 2-d theories this 
angle dependence can be predicted, but lattice animals 
are not conformally invariant. 

Apart from these specific problems we believe that the 
present simulations have demonstrated again the power 
of sophisticated sequential sampling methods, and of 
PERM in particular. Although there are certainly many 
problems where other MC strategies are more efficient, 
there are by now many examples where PERM seems un- 
challenged by any other known method. Unfortunately 
(or, rather, fortunately for the livelyhood of the subject) 
it is hard to predict when PERM or a similar strategy 
will be the method of choice. But we are confident that 
lattice animals will not be the last such problem. 
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